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(57) Abstract 



The present invention provides apparatus and methods for determining the concentration of sample components of a sample by an 
analytical technique that yields a spectrum that can be written as Y(w) - P(w) . C. The apparatus and methods of the invention account 
for experimental errors that give rise to distortions in the observed spectrum and mat consequently result in inaccurate determinations of 
sample component concentrations. The invention accounts for such errors by modeling the total experimental error as the sum of one or 
more types of errors that can be written as ( . K. The spectrum is then modeled asY-P.C + (.K. Using the observed spectrum, known 
values for P, and a mathematical model for ( , this equation can be solved for the best fit value of the sample component concentrations, C, 
and the magnitudes of the errors, K. The method can be used for any error that can be modeled in the foregoing manner, such as a shift in 
the spectrum. Particular types of shift include constant shift as well as linear shift across the entire spectrum. The apparatus and methods 
are advantageously used in absorbance spectroscopy and chromatography. 
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DETERMINATION OF COMPONENT CONCENTRATIONS TAKING ACCOUNT OF MEASUREMENT 
ERRORS. 



5 BACKGROUND 

Field of the Invention 

The present invention relates to the field of analytical spectrometry and, in 
particular, to apparatus and methods of correcting for determining analytical sample 

10 component concentrations that account for instrumental error. 
Summary of the Related Art 

A wide variety of analytical techniques have been developed over the years to 
detect and determine the concentrations of components of a sample. Some techniques are 
entirely spectrophotometry, such as ultraviolet- visible-infrared (UV-VIS-IR) absorbance 

15 spectroscopy and NMR spectroscopy. Other techniques, such as column chromatography, 
are not spectrophotometric per se f but often use spectrophotometry techniques to detect 
the presence of compounds. 

Each of these techniques generally provides a "spectrum" in which a dependent 
variable, typically the intensity of some quantity (e.g., absorbance), is plotted against a 

20 dependent variable (e.g., wavelength). Relative concentrations of sample components are 
determined by obtaining the best fit to the experimental spectrum by varying the relative 
spectral contribution of each component. This requires knowledge of the spectral features 
of the individual components. Owing to instrumental errors and/or deviations in 
experimental conditions, small shifts of the independent variable often occur that can lead 

25 to large changes in the relative concentrations of sample components determined from the 
measured spectra. To date, the prior art is devoid of suitable methods for determining the 
shift and correcting for it. Accordingly, new methods of accounting for shifts of the 
independent variable in analytical spectrophotometric techniques are desirable. 
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SUMMARY OF THE INVENTION 
The present invention provides apparatus and methods for the accurate 
determination of sample component concentrations. The apparatus and methods of the 
invention advantageously correct for experimental errors (including instrument induced 
5 errors) that would otherwise introduce errors into the measured sample component 
concentrations. The present invention compensates for a wide variety of phenomena 
exogenous to the analytical sample that contribute to and are manifested in the observed 
analytical spectrum, from which sample component concentrations are determined. 

The present invention can be employed advantageously in analytical 

10 spectrophotometry techniques such as UV-VIS-IR spectroscopy and NMR spectroscopy. 
The present invention can also be used in analytical techniques that are not 
spectrophotometry per se, but which incorporate spectrophotometry detection and/or yield 
a "spectrum-like" graph (Le. f a graph that resembles a spectrum). An example of a 
technique that yields a "spectrum-like" graph is column chromatography, in which 

IS materials are detected by absorbance spectroscopy at one (or more) wavelength, yielding 
a "spectrum-like" graph that depicts the intensity of the eluate absorbance as a function of 
time. As used herein, the term "spectrum" encompasses traditional spectrophotometry 
spectrum wherein a spectral intensity is plotted against radiation frequency, wavelength, 
or wavenumber (or some equivalent thereof), as well as "spectrum-like" graphs produced 

20 in techniques such as column chromatography. 

Under ideal conditions, an analytical spectrum can be described by the equation 

Y(o>) » P(o>) • C , (i) 
were to is a functional parameter (e.g., frequency, wavelength, wavenumber, or time), Y 
is a vector whose elements are the spectral intensities, C is a vector whose elements are 

25 the concentrations of the sample components, and P is a matrix whose elements are a 
measure of the magnitude of the contribution of each sample component to the spectral 
intensity at each value of u>. The elements of P are known quantities and can be, for 
example, the UV-VIS-IR extinction coefficients of each sample component at each 
wavelength A.. In practice, when one desires to use such analytical spectra to determine 

30 the concentrations of sample components that give rise to the observed spectrum, Y^, 
equation (i) can be used to obtain a best fit estimate of C by using Y^, in place of Y. 
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Owing to instrumental and other experimental errors, however, equation (i) 
frequently does not describe the observed spectrum well. In that case, equation (i) can be 
modified to incorporate a term representing the experimental error: 

Y = P • C + dY , (ii) 
5 where dY is the experimental error-induced deviation of the observed spectrum from the 
ideal. In the present invention, dY is written as: 

dY = 5'K , (iii) 
where K is a scalar representing the magnitude of the error and I is a vector whose 
elements are the relative errors at each value of o. Equation (ii) is then written as 
10 Y = P-C+5-K. (iv) 

The apparatus and methods of the present invention can be used to correct for experimental 
errors whenever the observed spectrum can be written in the form of equation (iv), i.e., 
whenever the experimental error can be written as in equation (iii). Equation (iv) is 
readily solved for the best fit sample component concentrations, C, and the magnitude of 
15 the error, K, to the observed spectrum, Y^, (using Y,*, in place of Y in equation (iv)) as 
described more fully below. 

In one aspect of the invention, the apparatus and methods model the error, dY, as 
a shift of the entire spectrum by an amount do> . In this aspect of the invention: 

dY = — do = Y'dto . (v) 
d(x> 

20 In this aspect of the invention, the observed spectrum is estimated as: 

= P • C + Y'dQ . (vi) 
Depending on the model chosen for the shift, du can be a scalar (i.e. , the same for all o>) 
or a vector whose elements vary with o>. 

In one embodiment of this aspect of the invention, using an appropriate model for 
25 do, equation (vi) takes the form of equation (iv), which is then solved for the best fit 
values of C and K. 

In another embodiment of this aspect of the invention, the apparatus and methods 
account for shift due to experimental error by adjusting the entire spectrum (to yield an 
adjusted spectrum, Y^) using some weighted average dw of previously determined 
30 values of the shift, do>. The entire spectrum is then adjusted by du> : 

Y^ = Y^u + dU ). (vii) 
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When dco is small, Y' da> is a good estimate of the shift in the spectrum, which is 
then preferably calculated from: 

Yadj = Y** + Y' dZ> . (viii) 
On the other hand f when dco is large, it is preferable to account for the shift by using 
5 equation (vii) directly. In either case, in this embodiment of the invention, is used in 
either equation (i) or (vi) to obtain the best fit value of C. 

In another embodiment of this aspect of the invention, the form of Y in equation 
(i) is used in equation (v) to obtain an expression for dY in terms of the derivative to P 
with respect to co: 

10 dY = —do) - ^^do) = — -Cdco ^ 

8co 8co dto 9 K } 

and the estimated spectrum is: 

Y nt = P • C + P' • C dco . ( X ) 
Equation (x) reduces to the form of equation (iv) when dco is modeled appropriately and 
can then be solved for C. 
15 In another embodiment of this aspect of the invention, the apparatus and methods 

of the present invention account for shift adjusting P using some weighted average of 
previously determined values of dco . An adjusted value of P that accounts for the shift is 
then given by: 

P^ = P (co + dco ) , (xi) 
20 If dco is large, equation (xi) is preferably used directly. If dco is small, however, 
(8P/dco)dco is a good estimate of the shift-induced change in P, and P^ is then preferably 
obtained from: 

P«ij = P + ^-d<o . (xii) 
dco 

In either case, P^ is then used in either equation (i) or (x) to obtain the best fit value of 
25 C. 

In this aspect of the invention, the shift can be modeled in a number of ways. In 
one embodiment, the shift is modeled as being constant across the entire spectrum. In this 
model dco = S, a scalar. In another embodiment, the shift is modeled as varying linearly 
about a central value of co, g> c , and is given by the vector whose elements are 
30 dco i =(co i -co c )M , where M is the magnitude of the shift. In this embodiment, dco is a 
vector. When M>0 the entire spectrum is magnified about to c ; when M<0, the entire 
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spectrum is compressed about g> c . In yet another embodiment, the experimental shift is 
modeled as a linear combination of the two foregoing models. 

As described more fully below, each of these models can account for errors in 
observed spectrum that are manifested by shifts in o>. The apparatus and methods of the 
5 present invention incorporate the foregoing equations to compensate for experimental shift 
in observed spectra and thereby yield more accurate measurements of the sample 
component concentrations. 

In a particularly preferred embodiment of the present invention, the apparatus and 
methods of the invention employ UV-VIS-IR absorbance spectroscopy. Co-oximeters, 
10 which measure the relative concentrations of blood components, are an example of one 
embodiment. In this embodiment, Y is the absorbance spectrum A, P is the matrix of 
extinction coefficients, E, and the independent variable, o, is the wavelength, A. 

In another preferred embodiment, the apparatus and methods of the invention 
employ chromatographic means. In this embodiment, ca is the time at which the sample 
IS components elute, Y is the absorbance at the wavelength at which the components are 
detected, and P is a matrix of the relative absorbances of each sample component as a 
function of elution time. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The apparatus and methods of the present invention are useful for correcting for a 
20 wide variety of instrumental and other experimental errors in a wide variety of 
spectrophotometry and analytical techniques. In general, the apparatus and methods of 
the present invention can be used in conjunction with any spectrophotometry and/or 
analytical technique yielding a spectrum that can be described by the equation: 

Y(u) = P(o>) • C (1) 
25 where Y is a vector whose tt nT components Y s are the intensities of the spectrum at each 
of the w m" values of the independent variable to (o^), C is a vector whose u n" elements 
are the concentrations of the sample components that contribute to the measured response 
Y, and P is an "m x n" matrix whose elements ? t] relate the contribution of component C . 
to the intensity Y;. "m" is an integer and equal to the number of values at which Y is 
30 measured. u n" is also an integer and equal to the number of sample components 
contributing to Y. In UV-VIS-IR spectroscopy, for example, Y. is the sample absorbance 
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at wavelength o>i, and the P i} are the extinction coefficients of absorbant "j" at wavelength 
ci> ( . In column chromatography, as another example, Y { is the intensity of the absorbance 
at time co,, and the are the extinction coefficients of the absorbant u j" at the monitoring 
wavelength at elution time o>j under the particular elution conditions (e.g. t buffer identity 
5 and strength, pH, temperature, etc.). In practice, when one desires to use such an 
analytical spectrum to determine the best estimate of sample components concentrations, 
equation (1) can be used by inserting for Y and solving for the best fit value of C. 

As will become clearer below, it will be appreciated by those skilled in the art that 
the present apparatus and methods are ideally suited to analytical techniques in which the 

10 identity of the "n" sample components and the matrix elements of P are known and it is 
desired to determine the best estimate concentrations of the sample components from the 
observed (measured) spectrum Y^. 

The present apparatus and methods can be employed advantageously in analytical 
spectrophotometric techniques such as UV-VIS-IR spectroscopy and NMR spectroscopy. 

15 The present methods can also be used in analytical techniques that are not 
spectrophotometric perse, but which incorporate spectrophotometric detection and/or yield 
a "spectrum-like" graph (i.e., a graph that resembles a spectrum). An example of a 
technique that yields a "spectrum-like" graph is column chromatography, in which 
materials are detected by absorbance spectroscopy at one (or more) wavelength, yielding 

20 a 44 spectrum-like" graph that depicts the intensity of the eluate absorbance as a function of 
time. As used herein, the term "spectrum" encompasses traditional spectrophotometric 
spectrum wherein a spectral intensity is plotted against radiation frequency, wavelength, 
or wavenumber (or some equivalent thereof), as well as "spectrum-like" graphs produced 
in techniques such as column chromatography. 

25 It is frequently the case that errors in the observed spectrum are introduced due to 

less than ideal instrument performance and/or experimental technique. As used herein, the 
term "experimental error" means any error {e.g., arising from less than ideal instrument 
performance or suboptimal experimental technique) resulting in a deviation of the 
measured spectrum from the theoretical ideal. These errors translate into errors in the 

30 measured sample component concentrations. These errors in the observed spectrum can 
be accounted for by rewriting equation (1) as 
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= P • C + dY , (2) 
where dY is the experimental error-induced deviation in Y^ from the ideal. The present 
apparatus and methods corrects for such errors by incorporating a mathematical model for 
dY: 

5 dY « ? • K (3) 

where K is a vector of dimension "r" whose elements are the magnitudes of the errors for 
each of the u r" types of error modeled as contributing to the spectrum, and $ is an "m x 
r" matrix comprised of tt r w vectors each of whose "m" elements are the relative errors at 
each value of co for each type of error modeled as contributing to the spectrum, "r" is an 
10 integer and equal to the number of types of errors modeled as contributing to the spectrum. 
Where only one type of error is modeled r = 1 , and K is a scalar and ^ is a vector. 
Equation (3) is then written as: 

+ . (4 ) 

In its broadest aspect, the apparatus and methods of the present invention incorporate 

15 corrections for instrumental and/or other experimental induced errors in measured spectra 
whenever the spectra can be written in the form of equation (4). 

In this aspect, therefore, the invention comprises an improved apparatus and 
methods of determining the concentrations of sample components using an analytical 
technique that yields a spectrum that can be estimated as equation (1), the improvement 

20 comprising correcting for experimental error by modeling the experimental error as "r" 
types of errors given by the product 5 - K, where K is a vector whose "r* elements are 
the magnitudes of each of the M r" types of experimental errors and £ is an w m x r" matrix 
whose elements are the relative errors at each value of o> for each type of experimental 
error, adding the product £ • K to the estimated spectrum as in equation (4), and solving 

25 for the best fit values of C and K. 

As used in the present invention r * 1; preferably r £ 10; and most preferably 1 £ 
r s 3. n * 1 and preferably 1 * n s 20. min + r and preferably m is about twice n + 
r. 

Equation (4) is readily solved for the best fit values of C and K by least squares 
30 analysis. To do so, equation (4) is "collapsed" by defining an "m x (n+r)" augmented 
matrix, P 5 , which has the form: 
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10 



15 



20 



;•■ P i.n 5 U 5 



(5) 



and an augmented vector, C K , of dimension "n+r," which has the form: 



c K - 



(6) 



Equation (4) then becomes: 

Y ob% = P r C K . 
The least squares solution to equation (7) is given by: 



C K = Pr* 



L obi * 



(7) 



(8) 



where 



P 5 * = (PJ • P<)' 1 • P 5 T (9) 

is the least squares transformation matrix. E.g., Noble and Daniel, Applied Linear 

Algebra, pp. 57-65 (Prentice-Hall, Inc., N.J., 1977). C K is readily determined from 

equation (8) using standard algorithms. See, e.g., Press et al. t Numerical Recipes: The Art 

of Scientific Computing (Cambridge University Press, Cambridge 1986). The first u n n 

25 elements of C K are the best fit values of the w n" sample component concentrations and the 

remaining "r n elements are the best fit magnitudes of the errors. 

In another aspect of the invention, the apparatus and methods model the error, dY, 

as a shift of the entire spectrum by an amount do> . In this aspect of the invention: 

dY = — do = Y'do . (10) 
do) 

30 In one embodiment of this aspect of the invention, the spectrum is estimated by substituting 
equation (10) into equation (2) to obtain 



Yob. - P • C + Y'dto 



(ID 
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The term Y'dco must be capable of being written in the form of equation (3) and, 
therefore, equation (11) in the form of equation (4). As will be demonstrated below, this 
can be so even when dco itself is other than a scalar. . In that instance, do> is decomposed 
into a scalar K and a vector or matrix and the vector or matrix combined with Y r to form 
5 a product matrix 5 . Equation (1 1) is then solved in the same manner described above for 
equation (4). Specific examples of this method are presented in greater detail below. 

In another embodiment of this aspect of the invention, the shift due to experimental 
error is accounted for by adjusting the entire spectrum using some weighted average of the 
magnitude of the shift, dco , obtained, for example, from previously determined values 
10 of the shift, dco. The entire spectrum is corrected for experimental error by shifting it by 
an amount dco : 

v.* = Y obs (co + dZT ) . (12) 
When dco is large relative to the resolution of the measured spectrum, equation (12) is 
preferably used directly to obtain the adjusted, corrected spectrum, Y^. When dco is 
IS small relative to the measured spectrum, Y' dco is a good estimate of the shift in the 
spectrum, which is then preferably calculated from: 

Ya.j = + Y' dco . (13) 
In either case, however, Y^ is then used in either equation (1), (11), or (19) (see infra) 
in place of and the equations solved to obtain the best fit value for C that is corrected 
20 for shift. 

dco is any suitable scalar representing the shift in co. In a preferred embodiment, dco 
is the previously calculated value of dco, or an average over the last 44 k" measurements, 
where "k" is 2 or more. In a preferred embodiment, K is 5 or more. In another preferred 
embodiment, k = 8. This method weights each of the last "k" measurements equally. 

25 Alternatively, a filter can be used to give greater weight to the most recent values of dco. 
In this embodiment, each of the last w k" values of dco are weighted by a factor w w,. " dco 
is then represented by a vector, dco", where each of the "k n elements dco^ is a previously 
determined value of dco, such that dco,* is the most recently determined value of dco and 
dco k " is the oldest value, dco is then obtained from the equation: 

30 dco = w T • do>\ (14) 
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wherein w T is the transpose of the vector w whose tt k w values "w^ are chosen such that 
w, * w 2 * . . w k and: 

E w i = 1 - (15) 
w can be determined in any suitable manner. For equal weighting of the last "k" values 
of dco, each W| is 1/k. Alternatively, where it is desired to give greater weight to the most 
recently determined values of dco, a function such as: 

w i = Znl^T^ (16) 




10 or 

1 /r- 1 

. . (17) 
i=i a 

can be used, wherein w a w is a real number greater than 0. Preferably w a n is greater than 
1 . Other suitable weighting functions are well known to those skilled in the art. 

In another embodiment of this aspect of the invention, the form of Y in equation 
15 (1) is used in equation (10) to obtain an expression for dY in terms of the derivative to P 
with respect to <o : 

dY = — dco = i^do) = — Cdto = P'-CdG> Q8) 
da> 6o> 6a> * 

and the best fit value of C corrected for shift is obtained from: 

Yo* » P - C + P' • C do) . (19) 
20 Equation (19) reduces to the form of equation (4) when dco is modeled appropriately and 
an estimated value for C, C« t is used in the expression P' * C. C est can be obtained, for 
example, from the solution to equation (1): 

C ttl = P 1 * Y^ (20) 

where 

25 P f = (P T • P)' 1 • P T (21) 

is the least squares transformation matrix and Y^ is the observed (or measured) spectrum. 
With this estimate equation (19) is then solved for the best fit value of C and the magnitude 
of the shift in the same way equation (4) is solved, as described above. 

In another embodiment of this aspect of the invention, the shift is accounted for by 

30 adjusting P using some weighted average of previously determined values of do>. An 
adjusted value of P that accounts for the shift is then given by: 
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= P(co + dco ) . (22) 
If dco is large, equation (21) is preferably used directly. If do is small, however, 
P' dco is a good estimate of the shift-induced change in P, and P^ is then preferably 
obtained from: 

5 = P + P' dco . (23) 

In either case, P^ is then used in either equation (1) or (19) in place of P to obtain a value 
for C that is corrected for shift. In either case, P^ is used in any of equations (1), (1 1), 
and (19) to determine the magnitude of the shift, dco, and a corrected value of C. 

The shift can be modeled in a variety of ways. For example, in one embodiment, 
10 the shift in co is modeled as being constant across the entire spectrum. In this embodiment 
dco is defined as a scalar S, which can be positive or negative, dco is given by some 
weighted average, S , of S. In this embodiment equation (11) takes the form: 

Y*, = P • C + Y' • S . (24) 
In this embodiment, the shift can be accounted for to obtain more accurate values of C by 
15 solving equation (24) in the manner described for solving equation (4), supra, where i = 
Y' and K=S. Alternatively, the entire spectrum may be shifted by amount S, preferably 
as in equation (12) when S is large: 

Y adj = Y*, (co + S) . (25) 
or as in equation (13) when S is small: 
20 Y^ = Y^ + Y'S. (26) 

In either case, however, Y adj is then used in either equation (1), (24), or (27) (see infra) 
in place of Y^ to obtain S and a best fit value of C that is corrected for the shift. 

Alternatively, the shift is accounted for using the derivative of the matrix P, as in 
equation (19): 

25 Y*. = P • C + P' • C S . (27) 

Equation (27) is solved in the same manner as equation (19) to obtain the magnitude of the 
shift S and a value of C that is corrected for the shift. Or, if a reasonable estimate of S 
is available, equation (22) can be used to correct for the shift by adjusting the matrix P, 

using: 

30 P gdj = P (co + S) . (28) 

if S is large, or: 

11 
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P«jj = P + P'S . (29) 
if S is small. In either case, P^ is then used in any of equations (1), (24), and (27) in 
place of P to obtain a value for C that is corrected for shift. If either equation (24) or (27) 
is used, a new value for the magnitude of the shift, S, is also obtained. 
5 In another embodiment for modeling experimental error arising from a shift in the 

spectrum, the shift is modeled as a compression/magnification of the spectrum about a 
central value of g>, g> c . The change in to due to magnification or compression in this 
model is given by 

da>i = (Qj - (D c ) M , (30) 
10 where is the "i*" component of the vector o>, whose elements are the values of o> at 
which measurements are taken, and M is the magnification/compression factor. It is seen 
that in this model do> is a vector. If M > 0, then the scale of the independent variable 0) 
is magnified. If M <0, then the scale is compressed. It is seen from equation (30) that 
the change in o> for which this method compensates is directly proportional not only to the 
IS magnification/compression factor M, but also to the distance from the central value of to, 
o) c . Thus, the greater the distance from the central value, Ci> c , the greater the change. 
In this embodiment equation (11) takes the form: 

Yob, = P C + Y' • A M . (31) 
where A is a diagonal matrix whose diagonal elements A H are (coj - co c ). Equation (31) 
20 is the same as equation (4), with I = Y' *A and K=M, and can be solved in the same 
manner. 

Alternatively, if a reasonable value M is available (as described in equations (14)- 
(17) and associated text), the spectrum may be shifted by amount A M , preferably as in 
equation (12) when M is large: 
25 Y^i = Y*, (o, + (u, - o> c )M ) , (32) 

or as in equation (13) when M is small: 

Y.-Y^ + Y'-AM . (33) 
In either case, however, Y adJ is then used in either equation (1), (31), or (34) (see infra) 
in place of Y obt to obtain M and a best fit value of C that is corrected for the 
30 magnification/compression type shift. 
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Alternatively, the magnification/compression type shift is accounted for using the 
derivative of the matrix P, as in equation (19): 

Yob, = P • C + A • P' ■ C M . (34) 
Equation (34) is solved in the same manner as equation (19) (by estimating C in the 
5 expression A • P' • C, which is equivalent to 5 in equation (4)) to obtain M and a value 
of C that is corrected for the magnification/compression type shift. Or, if a reasonable 
estimate of M is available (as described in equations (14)-(17) and associated text), 
equation (34) can be used to correct for the magnification/compression type shift by 
adjusting the matrix P, using: 
10 P^«P<»i + <tt|-«e)M). 05) 

if M is large, or: 

P^P+A P' C (36) 

if M is small. In either case, P^ is then used in any of equations (1), (31), or (34) in 

place of P to obtain a best fit value for C that is corrected for the 
15 magnification/compression type shift. If either equation (31) or (34) is used, a new value 

of M is also obtained 

In yet another aspect of the invention, the apparatus and methods incorporate the 

foregoing techniques to compensate both for shift and magnification/compression. In this 

aspect of the invention, the shift, do>, is given by 
20 dwj = S +(<*>, - o> c )M . (37) 

Using this expression for dco in equation (11) yields: 

Yo,,, = P • C + Y'S + A • Y' M . (38) 

Equation (36) can be solved in the very same manner as equation (4) by defining £ as the 

matrix [Y', A • Y'] and K as a vector [S,M]. Whereas the previous two embodiments 
25 illustrated models to equation (4) in which £ and K were a vector and a scalar, 

respectively, this embodiment illustrates a situation in which I and K are a matrix and a 

vector, respectively. 

In an alternative embodiment of this aspect of the invention, a weighted average of 
the scalars S and M are used to precalculate an adjusted spectrum, Y^, using equation (12) 
30 wherein da) is S + - g> c )M: 

= +S + (u>i - to c ) M) (39) 
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which is used if the sum S + (to; - coj M is large, or: 

» Y + Y' S + A • Y'M (40) 
when the sum S + (coj - M is small. In either case, Y^ is used in any of equations 
(1), (38) or (41) in place of Y^ to obtain new values for S and M and a best fit value for 
5 C that is corrected for both types of shift error. 

In another embodiment of this aspect of the invention, the change in ci> due to the 
combined effects of shift and magnification/compression can be accounted for by 
calculating the derivative of the matrix P. Using the results previously reported, equations 
(27) and (34), the estimated spectrum is: 
10 Y^ = P • C +P' • C S + A • P' ■ C M . (41) 

Equation (41) is analogous to equation (38) and can be solved in the same manner as 
equation (4) by defining the matrix I as [P' • C,A • P' • C] and the vector K as [S,M]. 

Alternatively and/or in addition, an adjusted matrix P v P^, can be calculated using 
previously determined values of S and M . Using equation (22), the adjusted matrix P is 
15 given by: 

P*ij,i = P(Wi +S + (o r o) c )M). (42) 
when the sum S + (o> - cj c ) M is large and by equation (23): 

= p + P' . C + A • P' • C. (43) 
when S -f (a) - o)^ M is small. In either case, P adj can be used in place of P in any of 
20 equations (1), (38), or (41) to obtain a best fit value of C that is corrected for both types 
of shift. If either equation (38) or (41) are used, new values for S and M are also 
obtained. 

In a particularly preferred embodiment of the present invention, the foregoing 
methods are applied to analytical absorbance spectroscopy. In this embodiment of the 
25 invention, equation (1) is the well known Beer-Lambert law: 

AeJX) = E(X) . C, (44) 
where A obs (X) is the absorbance spectrum measured as a function of the wavelength k, 
E(k) is the matrix of wavelength-dependent extinction coefficients for the absorbents, and 
C is the concentration of the absorbents. A is a vector whose w nT elements A 4 are the 
30 absorbencies at a m n discrete wavelengths A. if E is an M m x n w matrix whose elements E i} 



14 



WO 97/01751 PCMB96/00616 

are the extinction coefficients of component "j" at wavelength A.,, and C is a vector, each 
of whose "n" elements C, is the concentration of absorbant "j". 

Hence, it is seen that in this particularly preferred embodiment, the following 
correspondences exist: 

5 co = A. (45a) 

Y ot» = Kb, (45b) 

Y ' - t (45C> 
P = E . (45d) 

P' = ^ (45e) 

OA 

10. A| = (A, - k c ) (450 

These definitions are then used in equations (1) through (43) to measure the magnitude of 
the experimental error and calculate a value of C that is corrected for error, as described 
above. 

This embodiment is particularly useful, for example, in analytical UV-VIS-IR 
15 absorbance spectrophotometry (e.g. , as conducted with the Ciba-Conxing Diagnostics 800 
Series Co-oximeters) used to determine the concentrations of blood components in blood 
samples. The major components of blood that will ordinarily be used in the present 
methods are reduced hemoglobin (HHb), oxyhemoglobin (0 2 Hb), carboxy hemoglobin 
(COHb), methemoglobin (MetHb), sulfhemoglobin (SHb), and lipid. The wavelength- 
20 dependent extinction coefficients for hemoglobin-based blood components for use in E can 
be measured directly (preferably on a highly calibrated spectrophotometer) or obtained 
from the literature. E.g., Zijlstra et al., Clin. Chem. 37(9), 1633-1638 (1991). The lipid 
spectrum can be measured from intravenous fat emulsion (e.g., the commercially available 
lipid product intralipid) in an aqueous dispersion of about 10% by weight. In a 
25 particularly preferred aspect of the invention, the inventive method is used to correct for 
instrumental wavelength drift observed in UV-VIS-IR absorbance spectroscopy. 

In another embodiment, the foregoing methods are employed to correct for 
experimental error in column chromatography. In this embodiment, the independent 
variable a> is equal to the elution time, t. The spectrum is the magnitude of the signal 
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for detecting the presence of sample components in the eluate (e.g., absorbance, refractive 
index). The shape of the elution profile measured under standard and controlled conditions 
provides the elements of P. In this embodiment, a shift in elution time, dt, is equivalent 
to the shift in the independent variable, dco . Such a shift may arise due to deviations from 
5 the standard conditions of such parameters as the flow rate, solvent strength, and column 
temperature. The change in elution time can be modeled as a combination of scalar shift 
and linear shift, as described in equations (37M43) and associated text. For example, a 
change in column flow rate affects the elution time for unretained components and will 
change the elution time for retained components in approximate proportion to their initial 

10 elution times less the elution time for the unretained components. Therefore, all peaks will 
be shifted by a fixed delay plus one proportional to the elution time difference. 

All of the foregoing mathematical manipulations can be conducted using standard 
software packages, such as MATHCAD (MathSoft, Cambridge, MA). 

In all of the methods of the present invention, a sample spectrum is generated from 

15 which the sample component concentrations can be determined. The sample component 
concentrations are then determined by employing one or more of the methods described 
previously. 

The apparatus according to the invention incorporate the methods according to the 
invention to more accurately determine the concentrations of sample components. 

20 Accordingly, the apparatus according to the invention comprises a means for generating 
a spectrum (as defined hereinabove) of an analytical sample, a means for detecting the 
spectrum, a means for recording the spectrum, and a means for manipulating the spectrum 
according to any of the methods described herein. Myriad means for generating, 
detecting, and recording a spectrum exist and are well known to those skilled in the art. 

25 E.g., Hobart H. Willard et al., Instrumental Methods of Analysis (7th Ed., Wadsworth 
Pub. Co., Belmont, CA, 1988). The means for manipulating the spectrum comprises any 
computer means that can run software embodying one or more of the foregoing methods 
for determining sample component concentrations corrected for experimental error from 
the measured spectrum. Of course, the practitioner will appreciate that the apparatus need 

30 not be an integrated unit. In a preferred embodiment, the apparatus is a spectrophotometer 
capable of generating a sample spectrum in the UV, VIS, or IR regions of the 
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electromagnetic spectrum. In another preferred embodiment, the apparatus comprises 
column chromatography equipment. 

The following examples are provided to illustrate certain embodiments of the 
invention and are not intended, nor should they be construed, to limit the invention in any 
manner. 

EXAMPLES 

All mathematical manipulations described herein were performed using the 
MATHCAD software by MathSoft (Cambridge, MA). Ciba-Corning Diagnostic 
proprietary extinction coefficient matrices were used. 

Example 1 

Estimating Wavelength Shift Using the Derivative of the Absorbance Spectrum 
The absorbance spectrum of eleven blood samples was measured at a resolution of 
1 nm and the fractions of HHb, 0 2 Hb, COHb and MelHb calculated using a least squares 
solution of equation (42). The results are as follows: 

0.5 -0.5 0.2 1.2 0.4 0 0.8 0.4 2.2 0.2 -0.4 

93.3 100.4 98.7 97.4 98.4 100 98 98.4 97.1 94.4 99.6 

5.1 -0.2 -0.1 0 0 -0.3 -0.1 0.2 0.6 4.8 0.4 

1.2 0.4 1.1 1.3 1.2 0.4 1.3 1 0.1 0.6 0.3 



(a) 



where each column is one sample and the rows are the fractional concentrations of HHb, 
0 2 Hb, COHb and MetHb, respectively. The data were shifted by 0.1 nm by fitting the 
measured spectrum with a cubic spline and shifting the fitted spectrum. The fractional 
concentrations were then determined by solving equation (42) for C using the least squares 
method. The result was: 

-0.3 -1.3 -0.5 0.4 -0.4 -0.8 0.1 -0.4 1.5 -0.5 -1.2 
93.7 101 99.3 98 99 100.6 98.6 98.9 97.6 94.9 100.2 
6.1 0.7 0.8 1 0.9 0.6 0.8 1.1 1.5 5.8 1.4 
0.5 -0.3 0.4 0.7 0.5 -0.4 0.6 0.3 -0.6 -0.1 -0.4 



(b) 
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The difference between the original data (a) and the shifted data (b) was determined to be: 
' -0.8 -0.8 -0.8 -0.8 -0.8 -0.8 -0.8 -0.8 -0.8 -0.8 -0.8 1 
0.4 0.6 0.6 0.5 0.6 0.6 0.6 0.6 0.5 0.4 0.6 
1 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 1 1 
0.7 -0.7 -0.7 -0.7 -0.7 -0.7 -0.7 -0.7 -0.7 -0.7 -0.8 



(c) 



As can be seen from matrix (c), the relative change in fractional concentrations due to a 
frequency shift as small as 0.1 nm can be quite pronounced, particularly for components 
that are present in relatively small concentrations. 

Using the derivatives of the absorption data spectra, the shifts and concentrations 
were calculated on a sample by sample basis using equation (22). The results for the shift 
were: 

(0.11 0.04 0.07 0.14 0.09 0.15 0.11 0.11 0.17 0.09 0.05), (d) 
and results for the concentrations were: 



0.6 -1 0 1.5 0.3 0.3 0.9 0.5 2.8 0.2 -0.7 

93.2 100.7 98.9 97.2 98.5 99.7 97.9 98.3 96.7 94.4 99.9 

5 0.3 0.2 -0.3 0.1 -0.8 -0.2 0.2 0 4.8 0.9 

1.2 0 0.9 1.6 1.2 0.7 1.4 1.1 0.6 0.5 -0.1 



(e) 



The difference between matrices (a) and (e) was calculated to be: 

20 f 0.1 -0.5 -0.2 0.3 -0.1 0.4 0.1 0.1 0.5 0 -0.4 

-0.1 0.3 0.2 -0.2 0 -0.3 -0.1 -0.1 -0.4 0 0.3 

-0.1 0.6 0.3 -0.3 0.1 -0.4 -0.1 -0.1 -0.6 0.1 0.5 

0.1 -0.4 -0.2 0.3 -0.1 0.4 0.1 0.1 0.5 0 -0.4 



(0 



By comparing matrices (c) and (0 « is seen that the differences from the original, unshifted 
25 spectrum are much smaller when one accounts for the shift. The average shift from matrix 
(d) is 0.104 nm. Note that even though the individual estimates of S can vary by as much 
as 50%, the average value is exceptionally good. 
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Using the average shift of 0.104 nm to compensate for the shift in the spectrum of 
each sample, the fractional concentrations were calculated using equation (24): 

0.5 -0.5 0.3 1.3 0.4 0 0.9 0.4 2.3 0.3 -0.3 

93.2 100.4 98.7 97.4 98.4 100 98 98.3 97.1 94.4 99.6 

5 -0.3 -0.1 0 -0.1 -0.4 -0.1 0.2 0.6 4.7 0.4 

1.2 0.4 1.1 1.4 1.3 0.4 1.3 1 0.1 0.6 0.3 



(g) 



and the difference between the original fractional concentrations (matrix (a)) and the 
fractional concentrations calculated using the average over 11 values of the shift, S, 
(matrix (h)), was determined to be: 



(h) 



0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


-0.1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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Given that the average value of the shift was a very good estimate of the actual shift, it is 
not surprising that the calculated estimates of the concentrations matched the concentrations 
of the unshifted data so well. 

Example 2 

Estimating Wavelength Shift Using the Derivative of the Extinction Coefficient Matrix 
Using the same data that resulted in matrix (a) in Example 1 and shifting the 
spectrum by 0.1 nm, we compensated for the shift by the method of calculating the 
derivative of the extinction coefficient matrix as described in equations (18)-(27) and 
associated text using estimated values of the component concentrations obtained by solving 
equation (42). The estimated fractional concentrations were determined to be: 
/ \ 



0.6 -1.0 0.0 1.5 0.3 0.4 0.9 0.5 2.7 0.2 -0.8 

93.2 100.8 98.9 97.6 98.4 99.7 97.9 98.4 96.7 94.4 99.9 

5.0 0.3 0.2 -0.3 0.1 -0.7 -0.2 0.1 0.0 4.9 0.8 

1.3 0.0 0.9 1.6 1.1 0.8 1.4 1.1 0.6 0.6 0.0 



0) 
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which differs from the original concentrations of matrix (a) as follows: 

0.1 -0.5 -0.2 0.3 -0.1 0.4 0.1 0.1 0.5 0 -0.4 

-0.1 0.4 0.2 -0.2 0 -0.3 -0.1 0 -0.4 0 0.3 

-0.1 0.5 0.3 -0.3 0.1 -0.4 -0.1 -0.1 -0.6 0.1 0.4 

0.1 -0.4 -0.2 0.3 -0.1 0.4 0.1 0.1 0.5 0 -0.3 



0) 



(1) 



20 



The wavelength shifts* S, were determined to be: 

(0.1 0.04 0.07 0.13 0.09 0.15 0.11 0.11 0.16 0.09 0.05 0.03) . (k) 
The average shift, S , is 0.097 nm. 

Using the average estimated shift, S , and solving for the concentrations C using 
equation (74) yielded the following concentrations: 

0.5 -0.6 0.2 1.3 0.4 -0.1 0.9 0.4 2.3 0.3 -0.4 

93.2 100.4 98.7 97.3 98.4 100 97.9 98.4 97 94.3 99.6 

5.1 -0.2 -0.1 0 0 -0.3 -0.2 0.3 0.6 4.8 0.5 

1.2 0.4 1.1 1.4 1.3 0.3 1.3 1 0.1 0.6 0.3 ; 

which resulted in the following differences from the original concentrations in matrix (a): 

000 000 0.1 0000 

-0.1 0 0 -0.1 0 0 -0.1 0 -0.1 -0.1 0 

0.1 0 0 0 0 0.1 -0.1 0 0 0 0 

0 0 0 0.1 0 0 0.1 0 0 0.1 0.00 



(m) 



Example 3 

Averaging S Over the Previous Eight Measurements 
Using the data from Example 1 and focusing on the ninth sample, one can average 
over the first eight samples to obtain a value of S . The original fractional values for the 



25 ninth sample (obtained from the unshifted spectrum) were: 



2.2 1 
97.1 
0.6 
0.1 , 



(O) 
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(P) 



and the fractional concentrations after shifting the spectrum by 0.1 nm were: 

' 1.5^ 
97.6 
1.5 

The average shift of the first eight samples was 0.096. The fractional concentrations and 
wavelength shift were calculated using the derivatives of the extinction coefficients. The 
fractional concentrations obtained were: 



97 
0.6 
1 0-1 



(q) 



and the calculated wavelength shift was 0.06. The shift estimated for the first sample was 
0.1. Continuing to average the shift over the last 8 samples, the new average shift S is 
given by 

0 096x8 - 0 1 ♦ Ml = 0.103 . 
7 8 

Using this revised figure for S , the adjusted value for the extinction coefficient matrix, 

Ertj, was calculated and used to obtain a revised set of fractional concentrations: 

'2.3^ 
97 
0.6 



(r) 



which, for this sample, differs insignificantly from the previous estimate. 
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We claim: 

1 . An improved method of determining the concentrations of one or more components 
of an analytical sample from an observed spectrum that can be estimated by 

Y^to) - P(o>) • C, 

5 wherein Y^isa vector whose w nT elements are the magnitudes of the observed spectrum 
at each value of an independent variable to f C is the vector whose "n" elements are the 
estimated concentrations of u n n components that contribute to the measured spectrum, and 
p is a a m x n" matrix whose elements are the magnitudes of the contribution to the 
spectrum of each of tt n" components at each of the tt nT values of the independent variable 

10 o>, wherein the method comprises generating a sample spectrum from which the 
concentrations of the sample components can be determined and determining the sample 
- component concentrations from the spectrum, the improvement comprising correcting for 
experimental error by modeling the experimental error as w r ff types of errors given by the 
product £ ■ K, where K is a vector whose w r w elements are the magnitudes of each of the 

15 w r w types of experimental errors and £ is an "m x r n matrix whose elements are the 
relative errors at each value of <*> for each type of experimental error, adding the product 
I • K to the estimated spectrum as 

YoJco) = P((o) • C + 5-K 
and solving for the best fit values of C and K, wherein M n w and M r" are integers each 

20 greater or equal to 1 and tt nT is an integer at least w n+ r\ 

2. The method according to claim 1, wherein the best fit values of C and K are 
determined by least squares analysis. 

3. The method according to claim 2 wherein r is 1 and the experimental error is 
modeled as a shift in o> by amount do, and 5 = Y' = (9Y obt )/(3a>). 

25 4. The method according to claim 2 wherein r is 1 and the experimental error is 
modeled as a shift in a> by an amount dco, £ = P' • C, and 

5 . An improved method of determining the concentrations of one or more components 
30 of an analytical sample whose observed spectrum can be estimated by the equation 

YoJu) = P(o>) * C 
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wherein Y^ is a vector whose u m n elements are the magnitudes of the observed spectrum 
at each value of an independent variable a>, C is the vector whose M n w elements are the 
estimated concentrations of w n" components that contribute to the measured spectrum, and 
P is a tt m x n" matrix whose elements are the magnitudes of the contribution to the 

5 spectrum of each of w n" components at each of the u m w values of the independent variable 
wherein the method comprises generating a sample spectrum from which the sample 
component concentrations can be determined and determining the sample component 
concentrations from the spectrum, the improvement comprising correcting for experimental 
error by modeling the experimental error as a shift of the spectrum by an amount do, 

10 estimating a shift, dco , calculating an adjusted spectrum, Y^, using either the equation: 

Y.4 - Y^co 4- d^ ) 



20 



30 



or: 



Yadj = Y^ + Y' • do 



and determining the concentrations, C, by solving either the equation: 
15 = P • C 

or: 

Y^ = P • C + I • K. 
for the best fit value of C, where I takes the form £ = Y' = (9Y obs )/(aa>) or 5 = P' ' C, 
where 



d(0 



u m" and u n n are intergers, tt n" is greater or equal 1 and "m" is at least "n." 
6. A method according to claim 5 wherein do> is the weighted average of u k n 
previous values of da>, wherein k is an integer greater or equal to 1, and do> is 
determined from the equation 
25 dco =w T -dco*, 

where w T is the transpose of a vector w of length u k" whose elements are the relative 
weights to be applied to each of the previous tt k w values of do, do* is a vector of length 
M k w whose elements are previously determined values of dco, and w satisfies the equation: 

i>,«i • 



i'l 



The method according to claim 6, wherein each element of w is 1/k, or: 
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a*i/ i=i a*i 




or 



where "a" is any real number greater than 1 and M i w is an integer from 1 to u k" 

8. An improved method of determining the concentrations of one or more components 

of an analytical sample whose observed spectrum can be estimated by the equation 

10 YeJco) = P(g>) • C 

wherein is a vector whose "m" elements are the magnitudes of the observed spectrum 
at each value of an independent variable g>, C is the vector whose "n" elements are the 
estimated concentrations of "n" components that contribute to the measured spectrum, and 
P is a "m x n" matrix whose elements are the magnitudes of the contribution to the 

15 spectrum of each of "n" components at each of the "m" values of the independent variable 
(o, the improvement comprising correcting for experimental error by modeling the 
experimental error as a shift in o> by an amount dco, estimating a shift, do> , calculating 
an adjusted P matrix, P^, using either the equation: 

P^ = P obs ((o + dco ) 

20 or: 

P« Jj = Po ta + P'- . 
and determining the concentrations, C, by solving either the equation: 

Yobs = I\dj " C 

or: 

25 Y^(u>) = P idj (G>) • C + 5 ' K. 

for the best fit value of C, where £ takes the form I = Y' = (dY«Jf(du>) or I « P' C, 
where 

r - *■ 

and "n" is greater or equal to 1 and "m" is at least "n." 
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9. A method according to claim 8 wherein the weighted average of dco, dco t is the 
weighted average of k previous values of do>, wherein k is greater or equal to 1, and 
determined from the equation 



wherein w T is the transpose of a vector w of length "k" whose elements are the relative 
to be applied to each of the previous k values of dw, wherein M k n is an integer greater or 
equal to 1, da>* is a vector of length "k* whose elements are the previously determined 
values of d&> , and w satisfies the equation: 



10. The method according to claim 9, wherein each element of w is 1/k, or: 



wherein tt a" is any real number greater than 1, and M i" is an integer from 1 to w k." 

11. The method according to any of claims 3 t 4, 6, or 9, wherein the shift in co is 
modeled as constant across the entire spectrum and equal to a scalar .do = S. 

12. The method according to any of claims 3, 4, 6, or 9, wherein the shift in a> is 
modeled as dcoj = (g>; - co^M, wherein dcOj is the shift at the i* value of co, co, is the i* 
value of co, a> c is a constant value of co, and M, which can be any real number, is the 
magnitude of the error. 

13. The method according to claim 11, wherein Y is an absorbance spectrum of the 
sample, P is the extinction coefficient matrix for the absorbing components of the sample, 
and co is the frequency at which the absorbance is measured. 

14. The method according to claim 12, wherein Y is an absorbance spectrum of the 
sample, P is the extinction coefficient matrix for the absorbing components of the sample, 
and co is the frequency at which the absorbance is measured. 

15. An apparatus for determining the concentration of components in a sample 
comprising 

a means for generating a spectrum from the sample; 
a means for detecting the spectrum; 



dco = w T • dco* , 



Ev 



or 
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a means for recording the spectrum; and 

a means for determining the concentrations of sample components from the 
spectrum according to the method of claim 13. 

16. An apparatus for determining the concentration of components in a sample 
5 comprising 

a means for generating a spectrum from the sample; 
a means for detecting the spectrum; 
a means for recording the spectrum; and 

a means for determining the concentrations of sample components from the 
10 spectrum according to the method of claims 14. 

17. The apparatus according to claim 15, wherein the apparatus comprises a 
spectrophotometer capable of generating a sample spectrum in the UV, VIS, or JR regions 
of the electromagnetic spectrum. 

18. The apparatus according to claim 16, wherein the apparatus comprises a 
15 spectrophotometer capable of generating a sample spectrum in the UV, VIS, or IR regions 

of the electromagnetic spectrum. 
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